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Abstract 

Forest soils store vast amounts of terrestrial carbon, but we are still limited in 
mechanistic understanding on how soil organic carbon (SOC) stabilization or 
turnover is controlled by biotic and abiotic factors in forest ecosystems. We 
used phospholipid fatty acids (PLFAs) as biomarker to study soil microbial 
community structure and measured activities of five extracellular enzymes 
involved in the degradation of cellulose (i.e., P-l,4-glucosidase and cellobiohy- 
drolase), chitin (i.e., P-l,4-N-acetylglucosaminidase), and lignin (i.e., phenol 
oxidase and peroxidase) as indicators of soil microbial functioning in carbon 
transformation or turnover across varying biotic and abiotic conditions in a 
typical temperate forest ecosystem in central China. Redundancy analysis 
(RDA) was performed to determine the interrelationship between individual 
PFLAs and biotic and abiotic site factors as well as the linkage between soil 
microbial structure and function. Path analysis was further conducted to 
examine the controls of site factors on soil microbial community structure 
and the regulatory pathway of changes in SOC relating to microbial commu- 
nity structure and function. We found that soil microbial community struc- 
ture is strongly influenced by water, temperature, SOC, fine root mass, clay 
content, and C/N ratio in soils and that the relative abundance of Gram-nega- 
tive bacteria, saprophytic fungi, and actinomycetes explained most of the vari- 
ations in the specific activities of soil enzymes involved in SOC 
transformation or turnover. The abundance of soil bacterial communities is 
strongly linked with the extracellular enzymes involved in carbon transforma- 
tion, whereas the abundance of saprophytic fungi is associated with activities 
of extracellular enzymes driving carbon oxidation. Findings in this study dem- 
onstrate the complex interactions and linkage among plant traits, microenvi- 
ronment, and soil physiochemical properties in affecting SOC via microbial 
regulations. 



Introduction 

The interactions between above- and below-ground com- 
ponents play an important role in driving ecosystem pro- 
cesses, but their underlying mechanisms are yet poorly 
understood. Thus, there have been continued and growing 
interests to elucidate the explicit relationships between veg- 
etation and below-ground processes and to seek mechanis- 
tic understanding on the contribution of soil biota and 



associated processes to ecosystem functioning (Wardle 
et al. 2004; Xiao et al. 2007; De Deyn et al. 2008; Jin et al. 
2010; Zhou et al. 2013). 

Soil microorganisms are a critical link between shifts in 
the composition of dominant vegetation and fundamental 
shifts in ecosystem functioning (Waldrop et al. 2000; Pres- 
cott and Grayston 2013; Prescott and Vesterdal 2013). Soil 
active microbial communities play a central role in below- 
ground processes, especially in mediating soil organic 
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matter decomposition and nutrient cycling in forest ecosys- 
tems (van der Heijden et al. 2008). An improved knowl- 
edge on the interactions between site factors and soil 
microbial communities in facilitating forest soil organic 
matter decomposition and soil carbon sequestration is 
increasingly recognized as key to understanding feedbacks 
of terrestrial ecosystem processes to global climate change 
(Singh et al. 2010). However, owning to the complexity of 
below-ground processes as well as technical difficulties to 
experimentally manipulate soil microbial structures and 
activities, significant gaps remain in our current under- 
standing on how soil microbial communities are controlled 
by complex interactions of biotic and abiotic site factors 
and how the structural shifts in soil microbial communities 
are linked to alterations of their functioning, such as in 
mediating soil organic carbon (SOC) dynamics (Hackl 
et al. 2005; Brockett et al. 2012). While advancement in 
biotechnology has allowed for the development of tech- 
niques for structural analysis of soil microbial communi- 
ties, for example, the characterization of microbial 
functional components with phospholipid fatty acids 
(PLFA) analysis, determinations of the linkage between 
microbial structure and function still pose a critical issue 
(Prescott 2010). A lack of effective technical solution in 
directly relating microbial community structure to func- 
tioning in SOC transformation and turnover in natural 
ecosystems has been an obstacle for seeking mechanistic 
understanding on the feedbacks of changes in below- 
ground processes to ecosystem succession and functioning 
in response to climate and land-cover changes. 

The PLFAs extracted from soil are used as an indicator 
of microbial community structure, as certain groups of 
microorganisms have different "signature" fatty acids 
(Tunlid and White 1992). These lipids are only present in 
the membranes of viable organisms and are rapidly 
degraded in soil and as such are a measure of the organisms 
living at the time of sampling (Grayston and Prescott 
2005). A more advanced technique based on nucleic acid 
extraction and analysis has also been increasingly used for 
examining soil microbial community structure, particularly 
genes coding for ribosomal RNA (rRNA). However, the 
PLFA method is a rapid and inexpensive way of assaying 
the biomass and composition of microbial communities in 
soils and may even be more sensitive in detecting shifts in 
microbial community composition when compared to 
nucleic acid-based methods (Ramsey et al. 2006). Although 
the PLFA method cannot compete with the rRNA methods 
in the phylogenetic resolution by which a given community 
can be characterized, the latter provide little information 
on the phenotype and the activity of the microorganisms in 
the environment (Frostegard et al. 2011). 

Soil microbial community structure in natural ecosystems 
is known to be directly influenced by many edaphic factors, 



such as temperature, water, and soil physicochemical 
properties (Williams and Rice 2007; Collins et al. 2008; 
Angel et al. 2010; Castro et al. 2010). Vegetation type and 
structure, by modifying the site microclimate, the quantity 
and quality of litter, the production of root exudates, and/ or 
the allocation patterns of organic matter, directly and 
indirectly affect soil microbial community structure (De 
Deyn et al. 2008; Wardle et al. 2012; Jassey et al. 2013). Pre- 
vious studies have demonstrated that ecosystems may differ 
in the relative influences of site factors on soil microbial 
communities. For example, Brockett et al. (2012) found that 
soil water was the dominant factor influencing the potential 
function of the soil microbial community across forests of 
varying biogeoclimatic zones. In contrast, many studies show 
that SOC is closely related to soil microbial community 
structure and function under different types of vegetation 
(Grayston and Prescott 2005; Yao et al. 2006; Franklin and 
Mills 2009; Katsalirou et al. 2010). There are also findings 
that specific soil chemical properties, such as soil C/N ratio 
(Fierer et al. 2009), nutrient status (Lauber et al. 2008), and 
soil pH (Rousk et al. 2009), are highly correlated with soil 
microbial community composition and function, and that 
plant litter chemistry (Ushio et al. 2008; Strickland and 
Rousk 2010) and spatial pattern of soil properties (Ushio 
et al. 2010) can impose marked impacts on forest soil micro- 
bial function by changing soil microbial community compo- 
sition. Despite a much improved understanding on the 
controls of below-ground processes, however, there lacks a 
general model quantifying the relative contributions to, and 
the levels of control of, soil microbial processes when multi- 
ple biotic and abiotic site factors are involved. The underly- 
ing mechanisms on the functional linkages between plant 
communities and soil microbial communities and the 
impacts on forest biogeochemical processes that arise as a 
consequence of below-ground activities are much compli- 
cated by the multiple effects of various site factors, biotic or 
abiotic, leaving alone the complexity in interactions and 
autocorrelations among the factors that potentially lead to 
nonlinearity of biotic and environmental controls of ecosys- 
tem processes. Clearly, a holistic approach taking into con- 
sideration of all potential factors and drivers is necessary 
when examining the structure-function relationships of soil 
microbial communities in order to gain mechanistic under- 
standings on the controls and functioning of below-ground 
processes. 

In this study, using a typical temperate forest in central 
China as a model ecosystem, we examined the linkage 
between soil microbial structure and functioning in driv- 
ing SOC transformation and turnover and assessed the 
relative influences of biotic and abiotic site factors on soil 
microbial structure and function. Soil microbial commu- 
nity structure, as characterized by PLFA analysis, and 
microbial function in driving SOC transformation and 
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turnover, as represented by five extracellular enzymes 
involved in the degradation of cellulose (i.e., (3-1,4-gluco- 
sidase and cellobiohydrolase), chitin (i.e., P-l,4-N-acetyl- 
glucosaminidase), and lignin (i.e., phenol oxidase and 
peroxidase), were investigated across varying biotic and 
abiotic conditions established by setting up sampling plots 
in different types of forest stands naturally occurring or 
historically planted in the study area. Multivariate analysis 
and pathway analysis were used to identify the dominant 
biotic and abiotic site factors affecting soil microbial 
structure and function and to determine the regulatory 
pathway of changes in SOC relating to microbial commu- 
nity structure and function. Our objectives were to deter- 
mine that, in forest ecosystems, (1) how soil microbial 
communities are controlled by complex interactions of 
biotic and abiotic site factors; and (2) how and to what 
extend the structural variations in soil microbial commu- 
nities are linked to functioning in driving SOC transfor- 
mation and/or turnover. We anticipate that the outcome 
of the study would help with gaining better mechanistic 
understanding on the microbial regulation of SOC 
dynamics. We hypothesized that soil microbial commu- 
nity controls SOC stabilization or turnover through a 
linkage between structure and function and that the role 
of biotic and abiotic site factors in affecting SOC stock is 
operated via regulations on the structural attributes of soil 
microbial community (Fig. 1). 



Materials and Methods 
Study sites 

The study was located in field sites at the Baotianman 
Long-Term Forest Ecosystem Research Station in the Bao- 
tianman Nature Reserve (latitude 33°20'-33°36'N, longi- 
tude 111°46'-112°04'E, and elevation 600-1860 m a.s.L), 
in the east of the Qinling Mountains in central China. It 
is in a transitional zone from warm temperate to north- 
ern subtropical climatic region. The climate is of a conti- 
nental eastern monsoon type, with annual mean air 
temperature of 15.1°C and mean annual precipitation of 
900 mm (Luan et al. 2011). Precipitation occurs mainly 
in the summer months June through August (55-62%; 
Liu et al. 1998). The soils are of dystric cambisols (FAO- 
UNESCO soil classification system) developed on weath- 
ered arenites. 

The zonal vegetation at the study sites is typically a 
deciduous broadleaved forests dominated by, with 
changes in topography and position on slopes, Quercus 
aliena var. acuteserrata Maxim., Quercus glandulifera var. 
brevipetiolata Nakai., and Quercus variabilis Blume., 
respectively, in the canopy layer. Pinus armandii Franch. 
plantations were established in the region for timber pro- 
duction around 1956, but some stands developed into 
mixed forests due to natural regeneration of Q. aliena 



Biotic/abiotic interactions & regulation Microbial controls over soil carbon 

of soil microbial community transformation & turnover 




Linkage between microbial 
structure & function 



Figure 1. A conceptual model illustrating the interactions between site biotic and abiotic factors and regulations of soil microbial community, 
linkage between soil microbial structure and function, and microbial controls over soil carbon stabilization and turnover via effects on enzymatic 
activities. 
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and implementation of forest protection schemes. There 
are also some azonal tree species, such as Acer mono 
Maxim., Toxicodendron verniciflnum (Stokes) F. Barkley, 
Carpinus cordata BL, and Populus davidiana Dode, etc., 
occurring in low frequency in the local forests. The 
understory is composed of predominantly woody plants 
Lauraceae obtusiloba BL, Rubus corchorifolius Linn. F, Vitis 
amurensis Rupr., Carpinus turczaninowii Hance., Celastrus 
orbiculatus Thunb., and Clematis florida Thunb. in the 
shrub layer, with infrequent occurrence of Pteridium 
aquilinum (Linn.) Kuhn var. latiusculum (Desv.) Underw. 
ex Heller and highly sparse herbaceous plants Carex sider- 
osticta Hance., Carex lanceolata Boott, Rodgersia aesculifo- 
lia BataL, Duchesnea indica (Andr.) Focke, and 
Adenophora axilliflora Borb. in the field layer. 

Four forest types, each dominated by Q. aliena, 
Q. glandulifera, Q. variabilis, and mixture of P. armandii 
and Q. aliena, respectively, and an age sequence of Q. 
aliena forest (~40, ~80, and >160 years), were included in 
this study. Three 20 m x 20 m plots were set up on sep- 
arate sites for each of the forest types and the age 
sequence stands in April 2011. The stands of similar age 
(~80 years) were selected for forest types representing 
Q. aliena, Q. glandulifera, and Q. variabilis. The mixed 
Q. aliena and P. armandii forest stands were aged 
>50 years. Stand age was obtained from forest manage- 
ment records and increment core samples. All sites are 
under natural conditions and have not experienced appar- 
ent anthropogenic disturbance in recent history. 

Field survey was conducted on each plot to collect data 
on tree species composition and density, size of individual 
trees > 5 cm DBH, canopy openness, litter mass on for- 
est floor, fine root biomass, and soil physicochemical 
properties. 

Field measurements and sampling 

Soil temperatures at 10 cm depth (T soi i 10 ) were measured 
using HOBO Data Loggers (TidbiT v2 Temp, Onset 
Computer Corporation, Cape Cod, MA) at 30 min inter- 
vals. The average soil temperature for each plot was calcu- 
lated for 4 weeks prior to soil sampling, which 
represented the critical time period for microbial response 
to varying environmental conditions (Bell et al. 2009). 
Plant fine roots and litter samples were collected follow- 
ing soil sampling. Litter samples (including undecom- 
posed plant tissues, partially decomposed duff, and O e 
and O a horizons) were collected on each plot at 12 ran- 
dom locations with a 20 cm diameter ring frame and 
oven-dried at 65°C to constant weight. Twelve root core 
samples were also collected on each plot by soil layer for 
0-5 cm, 5-10 cm, and 10-20 cm using a 15-cm-inner- 
diameter soil-corer. They were washed clean using 0.4- 



mm mesh bags and separated into coarse (>5 mm), med- 
ium (2-5 mm), and fine (<2 mm) roots before being 
oven-dried at 65°C to constant weight. All litter and root 
samples were measured for total dry mass and analyzed 
for concentrations of C and N, and the acid soluble frac- 
tion and insoluble residue (i.e., the "Klason lignin") (Par- 
ton et al. 2007; Meier and Bowman 2010). 

Soil samples were collected in mid-August 2012, at 12 
systematically arranged locations on each plot at a distance 
of 5 m from the plot center in the directions of 0, 30, 60, 
90, 120, 150, 180, 210, 240, 270, 300, and 330°. They were 
collected at depths of 0-5 cm, 5-10 cm, and 10-20 cm 
using a 10-cm-inner-diameter soil-corer. The samples were 
mixed to obtain one composite sample for each of the soil 
layers on each plot. All soil samples were immediately 
stored on ice in insulated containers upon collection and 
sieved to pass a 2-mm mesh after returning to laboratory. 
Each processed soil sample was divided into two portions 
and stored in a refrigerator at — 20°C before being analyzed. 
One was freeze-dried for PLFA analysis and the other for 
enzyme assays and physicochemical analysis. 

Altogether, 216 forest floor litter samples, 648 fine root 
samples, and 54 composite soil samples were collected 
and analyzed in this study. 

Measurements of soil physicochemical 
properties 

The gravimetric soil water content (% SWC) was calcu- 
lated from the mass loss after drying the samples at 
105°C to a constant weight, for at least 48 h. We used a 
hydrometer method for analysis of soil particle size distri- 
bution (Sheldrick and Wang 1993). The soil pH was mea- 
sured by mixing the soil sample with deionized water at a 
1:2.5 ratio (w/v; Liu et al. 2012). The supernatants were 
measured using a pH meter (HI-9125, Hanna Instruments 
Inc, Woonsocket, RI). Soil organic carbon (SOC) content 
was analyzed by K 2 Cr 2 07-H 2 S04 calefaction method 
(Nelson and Sommers 1982). Soil total nitrogen (TN) 
content was analyzed using Kjeldahl digestion procedure 
(Gallaher et al. 1976). Soil ammonium nitrogen (NH 4 -N) 
and nitrate nitrogen (N0 3 -N) were determined from 10 g 
(dry mass) of soil using 50 mL 2M KC1 extraction proce- 
dure and analyzed colorimetrically by a continuous flow 
analyzer (SEAL AA3, Norderstedt, Germany). Microbial 
biomass carbon (MBC) was measured by fumigation- 
extraction method, using 0.5 mol/L K 2 S0 4 as extracting 
agent; a conversion factors of 2.64 were used to convert 
extracted C to MBC (Vance et al. 1987). The dissolved 
organic carbon (DOC) was determined by directly mea- 
suring carbon concentration in the 0.5 mol/L K 2 S0 4 
extract solution without fumigation treatment (Wu et al. 
2011). 
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Phospholipid fatty acids (PLFA) analysis 

Microbial community structure was assessed by analyzing 
the composition of extractable ester-linked PLFAs, using 
the method was described by Bossio and Scow (1998). 
Concentrations of individual PLFAs were calculated based 
on 19:0 internal standard concentrations. The indicator 
PLFAs were used for classification of microbial commu- 
nity types. Bacterial community was represented by 
PLFAs il4:0, 15:0, il5:0, al5:0, il6:0, 16:lco7c, 17:0, al7:0, 
il7:0, cyl7:0, 18:1oj7c, and cyl9:0 (Frostegard and Baath 
1996; Liu et al. 2012). Gram-positive bacteria are com- 
posed of PLFAs il4:0, il5:0, al5:0, il6:0, al7:0, and il7:0 
(Liu et al. 2012; Nie et al. 2013), Gram-negative bacteria 
16:lra7c, cyl7:0, 18:lco7c, and cyl9:0 (Zelles 1999; Liu 
et al. 2012), actinomycete 10Mel6:0, 10Mel7:0, lOMe 
18:0 (Frostegard and Baath 1996; Bell et al. 2009; Brockett 
et al. 2012), saprotrophic fungi 18:2<n69c (Nie et al. 
2013), and arbuscular mycorrhizal fungi (AMF) 16:lco5c 
(Olsson 1999; Brockett et al. 2012; Nie et al. 2013). Other 
PLFAs such as 14:0, 16:0, 16:1 20H, 16:lco9c, 17:lco8c, 
and 18:lco9c were also used for analysis of the microbial 
composition (Liu et al. 2012). 

Soil extracellular enzyme assays 

We assessed microbial function by measuring activities of 
five soil extracellular enzymes involved in degrading cellu- 
lose (|3-l,4-glucosidase and cellobiohydrolase), chitin (P- 
1,4-N-acetylglucosaminidase), and lignin (phenol oxidase 
and peroxidase), respectively. The activities of P-glucosi- 
dase (BG; EC:3.2.1.21), cellobiohydrolase (CBH; 
EC:3.2.1.91), and N-acetyl-P-glucosaminidase (NAG; 
EC:3.2.1.30) were determined by the conventional p-ni- 
trophenol assays (Parham and Deng 2000; Baldrian 2009). 

Phenol oxidase and peroxidase activities were deter- 
mined using 1-3,4-dihydroxyphenylalanine (L-DOPA) as 
substrate (Sinsabaugh et al. 1993; Li et al. 2010). For phe- 
nol oxidase, the reaction mixture composed of 2 mL 
5 mmol/L L-DOPA solution and soil slurry (1 g fresh soil 
with 1.5 mL 50 mmol/L sodium acetate buffer), and per- 
oxidase activity assays received 2 mL of 5 mmol/L-DOPA 
and soil slurry (1 g fresh soil with 1.5 mL 50 mmol/L 
sodium acetate buffer), plus 0.2 mL of 0.3% H 2 0 2 . All 
total enzyme activities were expressed as umol/g soil/h. 
The results of all enzymatic assays were expressed on a 
dry-weight basis. There were three replicates for each soil 
sample. 

The specific enzyme activity, which is defined as the 
activity of a given enzyme on per unit of SOC basis, was 
calculated for each of the five soil extracellular enzymes 
for comparison of their relative contributions to the 
dynamics of SOC decomposition (Baldrian et al. 2013). 



Soil microbial respiration measurements 

Soil microbial respiration (MR) was measured by determin- 
ing C0 2 evolution over a 12 -day incubation following the 
procedure described by Jin et al. (2010). Specifically, 20 g 
(dry weight equivalent) of each soil sample was adjusted to 
60% of water-holding capacity (WHC) and placed in a 250- 
mL Mason jar. Respired C0 2 was absorbed in 10 mL 
0.1 mol/L NaOH solution suspended inside the jar. All jars 
were sealed with rubber septum and incubated at 25°C and 
for 12 days in the dark. The NaOH solutions containing 
dissolved C0 2 were then each titrated with 0.05 mol/L HC1 
solution to determine the amount of C0 2 evolved at 3 -day 
intervals. After each reading, the jars were left open for half 
an hour to ventilate and then resealed. 

The microbial metabolic quotient (qC0 2 ) was calcu- 
lated by: 



Data analysis 

Based on biologic measurements and field survey and 
monitoring, a total of 21 variables were derived and used 
in the analysis of site factors on soil microbial community 
structures, including the slope and cross-section area at 
breast height (CAB) for capturing potential influences of 
topography and tree stands; soil physiochemical proper- 
ties consisting of r soiU0 , %SWC, pH, SOC, DOC, MBC, 
TN, NH 4 -N, N0 3 -N, C to N ratio in soil (C/N soil ), and 
clay content (%Clay); biotic factors consisting of FR, litter 
mass, litter lignin content (%lignini itter ), root lignin con- 
tent (%lignin root ), litter C to N ratio (C/Ni itter ), root C to 
N ratio (C/N root ), litter lignin to N ratio (lignin/N litter ), 
and root lignin to N ratio (lignin/N root ). 

Redundancy analysis (RDA) was performed on data of 
relative abundance of individual PLFAs, expressed as % 
mole of the total, and site biotic and abiotic variables for 
determining the interrelationships between soil microbial 
composition and environmental conditions (Andersen 
et al. 2010; Thorns et al. 2010). We also used RDA to 
determine the linkage between soil microbial community 
types and the specific activities of soil extracellular 
enzymes involved in SOC acquisition and oxidation. The 
forward selection procedure in RDA, based on Monte 
Carlo permutation with 499 iterations, was performed to 
determine the most significant discriminating variables on 
the composition of individual PLFAs and the specific soil 
extracellular enzyme activities, and the significant vari- 
ables (P < 0.05) were used in the final analyses. Signifi- 
cance tests for RDA were made using CANOCO software 
for Windows 4.5 (Biometris-Plant Research international, 
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Wageningen, The Netherlands). All data were log-trans- 
formed before being used in the analyses (Thorns et al. 
2010; Ushio et al. 2010; Wu et al. 2011). 

Path analysis was conducted to examine the controls of 
biotic and abiotic factors on soil microbial composition 
and the regulatory pathway of changes in SOC relating to 
microbial community structure and function. This was 
made by constructing structural equation models, known 
as the path models, based on the conceptual model 
shown in Fig. 1. Structural equation modeling (SEM) is 
an advanced multivariate statistical technique that allows 
for hypotheses testing of complex path-relation networks 
(Grace et al. 2007). Path models were established and 
separately tested based on hypothetical connections 
between the dominant biotic and abiotic site factors and 
the structural attributes of soil microbial communities, as 
well as the linkage of the soil microbial community types 
with the function as represented by activities of soil extra- 
cellular enzymes involved in SOC transformation (i.e., C 
acquisition) and turnover (i.e., C oxidation). The soil 
temperature, soil water content (%), soil clay content 
(%), SOC content (%), and soil C to N ratio (C/N sofl ) 
were retained to represent microclimate and soil physico- 
chemical properties, and the fine root biomass to repre- 
sent plant factor. In the path model depicting the 
hypothesis on regulatory pathway of SOC dynamics, the 
MBC was considered as an early indicator of changes in 
SOC (Zhang et al. 2013); the abundance of bacteria and 
saprophytic fungi was retained to represent the structural 
attributes of soil microbial community, which correspond 
to an increase in C availability and recalcitrant compound 
decomposition, respectively (Strickland and Rousk 2010; 
Barcenas-Moreno et al. 2011; Cusack et al. 2011), while 
the specific activities of the hydrolytic enzyme P-glucosi- 
dase and the oxidative enzyme phenol oxidase were 
retained to represent microbial function driving the C 
acquisition and the oxidation of complex C compounds, 
respectively (Cusack et al. 2011). 

The normality of data distribution was examined for 
heteroscedasticity, and all bivariate relationships were 
checked for signs of nonlinearities before the SEM analy- 
sis. These analyses were performed by using the maxi- 
mum-likelihood estimation procedures of Amos 17.0 
(SPSS Inc., Chicago, IL). The adequacy of model fitting 
was determined by several tests, including the % 2 test 
(P > 0.05, CMIN/df < 2), the goodness of fit (GFI) (val- 
ues >0.8 and <1), the root square mean error of approxi- 
mation (RMSEA) (values <0.05), and the comparative fit 
index (CFI) (values >0.9) (Grace et al. 2010; lassey et al. 
2013; Zhang et al. 2013). To obtain the most parsimoni- 
ous model, the nonsignificant indicators and pathways 
were eliminated in the final path model. For each set of 
analysis, R 2 values were obtained for each dependent 



matrix, representing the proportion of total variance 
explained by the model (Grace et al. 2010). 

Results 

Controls of site factors on soil microbial 
composition 

Among the 21 biotic and abiotic site variables, nine were 
found to relate significantly with the microbial composi- 
tion, including % SWC, SOC, r somo , % Clay, FR, C/N soib 
C/Ni itter , NH 4 -N, and NO3-N (Table 1). The first two 
components of RDA axes explained 73.5% of the variance 
in the relationship between soil microbial composition 
and the nine selective site factors. 

More than half of the variations in the microbial com- 
position are explainable jointly by the nine factors (sum 
of Lambda-B in Table 1). Forward selection of the nine 
factors in the RDA ordinations showed that the microbial 
composition was primarily influenced by %SWC, SOC, 
T so iiio> % clay content, FR, and C/N soi i (Table 1). Despite 
the relatively high marginal effect, NH 4 -N had negligible 
conditional effect on soil microbial composition. 

The RDA ordination biplot shows specific associations 
between the dominant site factors and individual PLFAs 
(Fig. 2). % SWC and % clay form positive associations pre- 

Table 1. Marginal and conditional effects of site factors on soil 
microbial composition obtained from the summary of forward selec- 
tion in Redundancy analysis (RDA). 



Variables 


Lambda- 
A 1 


Lambda - 
B 2 




f-ratio 4 


Soil water content 


0.22 


0.22 


0.002 


14.54 


(%SWC, w/w) 










Soil organic carbon (SOC) 


0.19 


0.14 


0.002 


11.71 


Soil temperature at 


0.19 


0.07 


0.002 


6.21 


10 cm depth (7" soN10 ) 










Soil clay content (%Clay) 


0.16 


0.05 


0.002 


4.55 


Fine root mass (FR) 


0.16 


0.04 


0.014 


3.49 


Soil carbon to nitrogen 


0.13 


0.02 


0.020 


2.77 


ratio (C/Nsoii) 










Nitrate nitrogen (N0 3 -N) 


0.17 


0.02 


0.146 


1.63 


Litter carbon to nitrogen 


0.12 


0.02 


0.136 


1.62 


ratio (C/Niitter) 










Ammonium nitrogen 


0.18 


0.00 


0.962 


0.32 



(NH 4 -N) 

Bold values indicate significant effects at P < 0.05. 

1 Describe marginal effects, which show the variance explained when 

the variable is used as the only factor. 

2 Describe conditional effects, which show the additional variance each 
variable explains when it is included in the model. 
3 The level of significance corresponding to Lambda-B when perform- 
ing Monte Carlo test at the 0.05 significance level. 
4 The Monte Carlo test statistics corresponding to Lambda-B at the 
0.05 significance level. 



638 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



Y. You et a/. 



Microbial Control of Forest Soil C Change 



dominantly with the Gram-negative bacteria 18:lco7c and 
the Gram-positive bacteria il5:0, and negative associations 
with the saprophytic fungi 18:2co6,9c, the Gram-positive 
bacteria al5:0, and the nonspecific PLFA 14:0; whereas 
^soiiio; C/N soi i, and C/N litter are associated with the same 
PLFAs but in reversed pattern. FR is positively associated 
with the nonspecific PLFAs 18:lco9c and 17:lco8c and nega- 
tively associated with the nonspecific PLFA 16:lco9c. SOC, 
NH 4 -N, and N0 3 -N are positively associated with the 
Gram-positive bacteria il4:0 and il5:0, the Gram-negative 
bacteria cyl7:0 and cyl9:0, the arbuscular mycorrhizal 
fungi 16:lco5c, and the actinomycetes 10Mel7:0, and nega- 
tively associated with the nonspecific PLFAs 16:lra9c and 
14:0. There are large numbers of PLFAs in various micro- 
bial communities not apparently associated with any of the 
dominant site factors (Fig. 2). 

Linkage between soil microbial community 
structure and function 

In the RDA ordination biplot, axis 1 and axis 2 explained 
84.7 and 11.9%, respectively, of the variance in the relation- 
ship between the specific activities of soil extracellular 
enzymes and the abundance in the five soil microbial com- 
munity types (Fig. 3). The five soil microbial types jointly 



explained 66% of variations in the specific activities of the 
five soil enzymes assessed (sum of Lambda-B in Table 2), 
with Gram-negative bacteria, saprophytic fungi, and arbus- 
cular mycorrhizal fungi appearing to be the most dominant 
controls (Table 2). The abundance in saprophytic fungi is 
positively associated primarily with the specific activities 
of phenol oxidase and peroxidase, and secondarily with 
NAGase, whereas Gram-negative bacteria and arbuscular 
mycorrhizal fungi positively correspond with the abun- 
dance of P-glucosidase and cellobiohydrolase. The abun- 
dance of actinomycetes is negatively associated with 
NAGase (Fig. 3). 

Regulatory pathway of soil microbial 
community structure 

The path model on the controls of soil microbial commu- 
nity structure by dominant site factors well passed all the 
statistical tests on adequacy = 2.8, P = 0.944, CMIN/ 
df = 0.355; GFI = 0.985; CFI = 1.000; RMSEA < 0.001) 
and explained 72.9 and 68.7%, respectively, of the vari- 
ance in the abundance of the bacterial and saprophytic 
fungal communities. Based on the path analysis, the 
abundance of soil bacterial community is directly and 
positively controlled by FR, %SWC, and %clay; whereas 




Figure 2. Redundancy analysis (RDA) ordination biplot of individual 
phospholipid fatty acids (PLFAs) and dominant site factors. 7~ SO ino, soil 
temperature at 10 cm depth; C/N| itt er, litter carbon to nitrogen ratio; 
C/Nsoii, soil carbon to nitrogen ratio; FR, fine root mass; SOC, soil 
organic carbon; NH 4 -N, ammonium nitrogen; N0 3 -N, nitrate nitrogen; 
%SWC, gravimetric soil water content (w/w); % Clay, soil clay 
content. 



Figure 3. Redundancy analysis (RDA) ordination biplot of the 
abundance in five soil microbial community types and the specific 
activities of five soil extracellular enzymes. Sap, saprophytic fungi; G + , 
Gram-positive bacteria; G~, Gram-negative bacteria; AMF, arbuscular 
mycorrhiza fungi; Actino, actinomycetes; BG, p-glucosidase; CBH, 
cellobiohydrolase; NAG, N-acetyl-p-glucosaminidase; PO, phenol 
oxidase; PER, peroxidase. 
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Table 2. Marginal and conditional effects of soil microbial community 
types on specific activities of soil extracellular enzymes involved in 
SOC transformation and turnover obtained from the summary of 
forward selection in Redundancy analysis (RDA). 



Variables 


Lambda- 
A 1 


Lambda - 
B 2 


P 3 


F-ratio 


Gram-negative bacteria (G~) 


0.44 


0.44 


0.002 


40.72 


Saprophytic fungi (Sap) 


0.40 


0.13 


0.002 


15.66 


Actinomycetes (Actino) 


0.11 


0.08 


0.002 


10.87 


Arbuscular mycorrhiza 


0.07 


0.00 


0.546 


0.71 


fungi (AMF) 










Gram-positive bacteria (G + ) 


0.06 


0.01 


0.552 


0.62 



Bold values indicate significant effects at P < 0.05. 

1 Describe marginal effects, which show the variance explained when 

the variable is used as the only factor. 

2 Describe conditional effects, which show the additional variance each 
variable explains when it is included in the model. 
3 The level of significance corresponding to Lambda-B when perform- 
ing Monte Carlo test at the 0.0S significance level. 
4 The Monte Carlo test statistics corresponding to Lambda-B at the 
0.05 significance level. 

the abundance of saprophytic fungi community is directly 
and positively controlled by C/N soi i and negatively by % 
SWC and %clay (Fig. 4). 

Apart from the direct controls, the path analysis also 
shows strong indirect effects of site factors by interacting 
with each other, resulting in cascade effects on the abun- 
dance of soil bacterial and saprophytic fungal communi- 
ties (Table 3). Taking into consideration of both direct 
and indirect effects, % SWC exerts strongest total positive 
control on the abundance of soil bacterial community, 
followed by % clay and FR in descending order, whereas 
the abundance of soil saprophytic fungal community is 
most strongly and negatively affected by % clay content, 



followed by % SWC (Table 3). r soi i 10 has a strong 
indirect positive effect on the abundance of soil sapro- 
phytic community. 

Regulatory pathway of soil carbon 
transformation and turnover 

The path model on the regulatory pathway of soil C- 
related variables also well passed all the statistical tests on 
adequacy (x 2 = 4.0, P = 0.776, CMIN/df = 0.577; 
GFI = 0.979; CFI = 1.000; RMSEA < 0.001) and 
explained 75.3, 23.8, and 48.9%, respectively, of the vari- 
ance in SOC, MBC, and qC0 2 . The path analysis points 
to direct and strong positive control of C acquisition 
activity by abundance in soil bacterial community, 
whereas the C oxidative activity is directly controlled by a 
positive effect of abundance in soil saprophytic fungal 
community and by negative effects of abundance in soil 
bacterial community (Fig. 5). SOC imposes strong posi- 
tive effects on the abundance of both bacterial and sapro- 
phytic fungal communities as well as MBC and is 
positively regulated by the C acquisition activity and neg- 
atively regulated by the C oxidative activity. Both the 
abundance of soil bacterial community and C oxidative 
activity imposes direct and strong positive effects on 
qC0 2 . 

When considering also the indirect effects, the C acqui- 
sition activity is most strongly and positively controlled 
by the abundance in soil bacterial community, with SOC 
exerting a strong indirect positive effect; the C oxidative 
activity is mostly strongly and negatively controlled by 
the abundance in soil bacterial community, followed by C 
acquisition activity, with SOC exerting a strong indirect 
positive effect (Table 4). SOC is most strongly regulated 
by the total positive effect of C acquisition activity, 



= 0.303 



R : = 0.433 



T sou 10 






- 






r R 2 = 0.45 


C/N i0i] 




FR 



~ >R 



R- = 0 732 



oSWC 



t 




oClav 



0.259 



= 0.729 



Bacterial abundance 



R- = 0.687 



Saprophytic fungal 
abundance 



Figure 4. Path model depicting the regulatory 
pathway of the controls of abundance in 
bacterial and saprophytic fungal communities 
by dominant site factors. Numbers on arrows 
are standardized direct path coefficients. R 2 
value represents the proportion of total 
variance explained for the specific dependent 
variable. Dash-line arrows indicate negative 
effects. r sol |i 0 , soil temperature at 10 cm 
depth; C/N soi |, soil carbon to nitrogen ratio; FR, 
fine root mass; %SWC, gravimetric soil water 
content (w/w); % Clay, soil clay content. 
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Table 3. Direct, 


indirect, and total 


path coefficients of the controls 


of dominant site 


factors on abundance 


of bacterial 


and saprophytic 


communities. 




















State variables 












Factors 


Effect 


S (7soino) 


S (%SWC) 


S (C/N soi ,) 


S(FR) 


S (Bacteria) 


S (Sap) 


F (% Clay) 


Direct 


—0.658 


0.488 


0 


0 


0.292 


—0.271 




Indirect 


0 


0.297 


—0.513 


0.253 


0.400 


-0.454 




Total 


—0.658 


0.785 


-0.513 


0.253 


0.692 


—0.725 


F (Tsoino) 


Direct 


0 


—0.452 


0.396 


0.649 


0 


0 




Indirect 


0 


0 


0.146 


-0.391 


-0.072 


0.325 




Total 


0 


—0.452 


0.542 


0.258 


-0.072 


0.325 


F (%SWC) 


Direct 


0 


0 


-0.321 


0.866 


0.383 


-0.409 




Indirect 


0 


0 


0 


0 


0.340 


—0.083 




Total 


0 


0 


-0.321 


0.866 


0.723 


-0.492 


F (C/N soil ) 


Direct 


0 


0 


0 


0 


0 


0.259 




Indirect 


0 


0 


0 


0 


0 


0 




Total 


0 


0 


0 


0 


0 


0.259 


F(FR) 


Direct 


0 


0 


0 


0 


0.392 


0 




Indirect 


0 


0 


0 


0 


0 


0 




Total 


0 


0 


0 


0 


0.392 


0 



% Clay, soil clay content; 7 soi | 10 , soil temperature at 10 cm depth; %SWC, soil water content; C/N soi |), soil carbon/nitrogen ratio; FR, fine root 
mass; Bacteria, abundance of bacterial community; Sap, abundance of saprophytic fungal community. 



R ; = 0.603 



C acquisition 




Figure 5. Path model depicting the regulatory pathway of the controls of soil organic carbon (SOC), microbial biomass carbon (MBC), and 
microbial metabolic quotient (qC0 2 ) by the structural attributes of soil microbial community and specific activities of soil extracellular enzymes 
involved in C acquisition and oxidation. Numbers on arrows are standardized direct path coefficients. R 2 value represents the proportion of 
total variance explained for the specific dependent variable. Dash-line arrows indicate negative effects. C acquisition is represented by the activity 
of p-glucosidase; C oxidation is represented by the activity of phenol oxidase. 



followed by strong indirect effect of the abundance in soil 
bacterial community; C oxidative activity and MBC 
impose strong total negative effects on SOC. MBC is posi- 
tively regulated by total effects of SOC and C acquisition 
activity and indirect effect of the abundance in soil bacte- 
rial community and negatively by C oxidative activity. 
qC0 2 is strongly and positively regulated by the total 
effects of C oxidative activity and indirect effect of the 
abundance in soil saprophytic fungal community, and 
negatively by indirect effects of C acquisition activity and 



SOC and the total effects of the abundance in soil bacte- 
rial community and MBC. 

Discussion 

A sound mechanistic understanding on the controls of 
biotic, climatic, and edaphic factors on soil carbon 
dynamics is a prerequisite for achieving a realistic predic- 
tion of changes in the regional soil carbon sink in 
response to changing environment. It is generally 
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Table 4. Direct, indirect, and total path coefficients of the controls of soil organic carbon (SOC), microbial biomass carbon (MBC), and microbial 
metabolic quotient (qC0 2 ) by abundance of soil bacterial and saprophytic fungal communities via impacts on specific activities of extracellular 
enzymes involved in carbon acquisition and oxidation. 



State variables 



Factors 


Effect 


S (Bacteria) 


S (Sap) 


^ (^-acquisition) 


5 (^oxidation) 


5 (SOC) 


S (MBC) 


qC0 2 


r (.DdCierlci; 


Direct 


u 


1 mn 
— I .uzu 


n r;77 
U. D / / 


— U.4-ZU 


n 
u 


n 
u 


n a c 1 

U.-40 I 




Indirect 


U.UDO 


n 1 n/i 

— U.I U4 


n TQ1 

u.zy i 


n 7^7 

— U. / j / 


n Q7£ 

U.O / O 


U.OOO 


— U.O / D 




i otai 


U. DUO 


1 1 ~)A 
— I.I £-*+ 


U.OUO 


— I.I// 


n Q7£ 

U.O / O 


n ££Q 

U.OOO 


— U.-4-Z4 


r \?>ap) 


Direct 


n 
u 


n 
u 


n 
u 


U./o^- 


n 
u 


n 
u 


u 




Indirect 


— U.UOD 


U.U I 4 


— U.UjC 


U.U^+Z 


n lie 
— U.I ID 


n n£7 

— U.UO/ 


U.ZZD 




i oxai 


n n££ 

— U.UOO 


U.U I 4 


— U.UJo 


U. 3/0 


n 111: 
— U.I ID 


— U.Uo/ 


U.ZZD 


F (^-acquisition) 


Direct 


0 


0 


0 


— 0.360 


0.749 


0 


0 




Indirect 


0.591 


-0.121 


0.341 


-0.271 


0.274 


0.780 


-0.538 




Total 


0.591 


-0.121 


0.341 


-0.631 


1.023 


0.780 


-0.538 


F (^-oxidation) 


Direct 


0 


0 


0 


0 


-0.329 


0 


0.651 




Indirect 


-0.233 


0.048 


-0.134 


0.146 


-0.074 


-0.308 


0.145 




Total 


-0.233 


0.048 


-0.134 


0.146 


-0.403 


-0.308 


0.796 


F (MBC) 


Direct 


0 


0 


0 


0 


-0.241 


0 


-0.503 




Indirect 


-0.171 


0.035 


-0.099 


0.107 


-0.055 


-0.226 


0.106 




Total 


-0.171 


0.035 


-0.099 


0.107 


-0.296 


-0.226 


-0.397 


F (SOC) 


Direct 


0.578 


0.471 


0 


0 


0 


0.763 


0 




Indirect 


0.100 


-0.616 


0.408 


-0.445 


0.226 


0.173 


-0.441 




Total 


0.708 


-0.145 


0.408 


-0.445 


0.226 


0.936 


-0.441 



Bacteria, abundance of bacterial community; Sap, abundance of saprophytic fungal community; C acqui s ition , C acquisition activity (represented by 
the activity of p-glucosidase); C oxic | at i 0n , C oxidative activity (represented by the activity of phenol oxidase); MBC, microbial biomass carbon; SOC, 
soil organic carbon; qC0 2 , microbial metabolic quotient. 



understood that SOC in forest ecosystems originates pre- 
dominantly from decomposing plant litter and that the 
transformation and turnover of plant-derived organic car- 
bon are mediated primarily by soil microbial communi- 
ties (Singh et al. 2010; Schmidt et al. 2011). As microbial 
communities in soils act to either stabilize or turning over 
plant-derived organic carbon (Liang et al. 2011; Miltner 
et al. 2011), it has long remained a challenge to directly 
link the soil microbial attributes with the capacity of soil 
carbon sequestration. Despite a profound progress over 
the past two decades in elucidating the microbial commu- 
nity structure as affected by site or environmental condi- 
tions in forest soils (Hogberg et al. 2007; Ushio et al. 
2008; Thorns et al. 2010), associating the soil microbial 
structure with functioning in driving soil carbon transfor- 
mation and turnover still occurs as an unresolved issue. 
In this study, by measuring a suite of complex site factors 
and soil microbial-related attributes under field condi- 
tions and applying multivariate analysis techniques, we 
examined ecosystem-level controls on soil microbial com- 
munity structure and the linkage between microbial struc- 
ture and function in soil carbon processes. Results show 
that at ecosystem level, factors relating to soil microcli- 
mate and fertility are the major controls on soil microbial 
structure and that the soil bacterial community and sap- 
rophytic community play distinct roles in mediating soil 



carbon processes in temperate forest ecosystem of the 
region. 

Research to date has uncovered some of the regulations 
of soil microbial community composition by individual 
biotic and abiotic factors and/or their interactive effects 
(Bezemer et al. 2006; Thorns et al. 2010; Brockett et al. 
2012). However, results are inconsistent concerning the 
critical drivers of differences in soil microbial structure. 
While the study of Hogberg et al. (2007) shows that soil 
C to N ratio and pH are the major factors associated with 
variations in soil microbial community structure, the 
study of Mitchell et al. (2010) suggests that plant com- 
munity composition better predicts changes in microbial 
community composition than soil properties. Thorns 
et al. (2010) found that plant community directly and 
indirectly affected soil microbial community composition 
through influences on the amount and quality of litter 
and soil physicochemical properties. Here, we show that 
soil water, SOC, soil temperature, soil clay content, fine 
root mass, and soil C to N ratio are all significant drivers 
of variations in soil microbial community composition. 
We found that soil water was positively associated with 
the abundance in bacterial community, especially the 
Gram-negative bacteria community, whereas soil temper- 
ature was positively associated with the abundance in 
saprophytic fungal community and negatively with the 
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abundance in bacterial community. Both soil temperature 
and soil water have been identified as primary drivers of 
variations in soil microbial community structure in other 
studies (Hackl et al. 2005; Djukic et al. 2010; Brockett 
et al. 2012). 

Factors relating to soil fertility, that is, SOC, NH 4 -N, 
NO3-N, and soil C to N ratio, were also found to corre- 
late with soil microbial community structure. Bacteria are 
typically found in soils with rich SOC, whereas the rela- 
tive abundance of saprophytic fungal community tends to 
increase with decreasing soil fertility (Grayston et al. 
2004; Grayston and Prescott 2005; Franklin and Mills 
2009; Katsalirou et al. 2010; Wu et al. 2011). The C to N 
ratio in both soil and litter was strongly and positively 
associated with saprophytic fungi in this study, consistent 
with the findings of positive correlation between the soil 
C to N ratio and the fungal to bacterial ratio by others 
(e.g., Hogberg et al. 2007; Fierer et al. 2009). Our finding 
that high soil clay content favors the abundance of bacte- 
rial community and constrains the abundance of sapro- 
phytic fungal community is well supported by literature 
(e.g., Hogberg et al. 2007; Lamarche et al. 2007; Fierer 
et al. 2009). Fine root mass, as a biotic factor, was found 
to associate mainly with nonspecific PLFAs. The influence 
of plant fine root mass on soil microbial community 
structure has been considered to involve several mecha- 
nisms, including alteration of soil microclimate, for exam- 
ple, changing soil water by uptake activity or modifying 
soil physical structure by root penetration, etc. (Prescott 
and Grayston 2013), and production of root litter (Lo- 
ranger-Merciris et al. 2006; Ushio et al. 2008) and root 
extrudes (Kuzyakov and Donmanski 2000; Jones et al. 
2004) as energy sources of soil microbial communities. 

Several recent studies have demonstrated that soil pH 
is among the primary site factors determining soil micro- 
bial community structure (Hogberg et al. 2007; Sinsab- 
augh et al. 2008). But in the current study, soil pH was 
not identified as a significant explanatory variable in the 
variations of soil microbial community composition. The 
narrow range of soil pH at our study sites (i.e., pH value 
4.14-4.73) and inclusion of more dominant site factors 
may contribute to this discrepancy. 

The nine dominant site factors explained only 58% of 
the total variations in soil microbial community composi- 
tion. The relatively low predictive power overall suggests 
there are other unaccounted factors influencing the 
microbial community composition. The unexplained 
variances may possibly due to the interference of root 
exudates and growth regulators (Pires et al. 2012), con- 
densed tannins (Ushio et al. 2010), total phenolic 
compounds (Ushio et al. 2010), and factors associated 
with soil type and structure (Girvan et al. 2003; Wakelin 
et al. 2008). 



Changes in the soil microbial community composition 
are likely to result in shift in the function of microbial 
communities, thus altering soil biologic and physiochemi- 
cal processes (Waldrop and Firestone 2006). Increases in 
the abundance of soil microbial communities that 
produce hydrolytic enzymes facilitate C acquisition in 
support of soil microbial primary metabolism (Cusack 
et al. 2011), whereas increases in the oxidative enzyme- 
producing microbial communities, mainly as saprophytic 
fungi, help with degrading complex compounds (Sylvia 
et al. 2004). An explicit understanding on the relationship 
between microbial community structure and function is 
crucial for predicting how shifts in microbial community 
structure may lead to changes in soil processes (Weand 
et al. 2010). In particular, identification of specific soil 
microbial attributes driving SOC transformation, and 
turnover is a significant step forward for gaining mecha- 
nistic understanding on the control of soil carbon sink 
and carbon sequestration potential in forest ecosystems 
(Lucas et al. 2007; Acosta-Martfnez et al. 2010). With the 
help of RDA ordination, which is a well-known multivari- 
ate analysis technique for discriminating the relative 
effects when multiple factors are involved in regulating a 
specific process, we were able to establish the relative 
linkage between soil microbial community types and 
functional-specific soil extracellular enzymes. We found 
that the relative abundance in the Gram-negative bacteria 
and arbuscular mycorrhizal fungi communities was asso- 
ciated with the specific activities of hydrolytic enzymes 
involved in C acquisition (i.e., (3-glucosidase and cellobio- 
hydrolase), whereas the relative abundance in the sapro- 
phytic fungi community was correlated with the specific 
activities of soil enzymes involved in degradation of chitin 
(i.e., NAGase) and lignin (i.e., phenol oxidase and peroxi- 
dase). Gram-negative bacteria generally grow fast to com- 
pete for simple organic compounds (Waldrop et al. 2000) 
and were found in this study to significantly relate to the 
specific activity of P-glucosidase involved in the break- 
down of complex organic compounds into small mole- 
cules substrates in favor of C acquisition through 
microbial community growth. Gram-positive bacteria 
were positively correlated with cellobiohydrolase responsi- 
ble for degrading complex compounds (Waldrop et al. 
2000; Bell et al. 2009). The saprophytic fungi possess the 
physiological capacity to produce phenol oxidase and per- 
oxidase essential in lignin depolymerization (Colpaert and 
Laere 1996; Courty et al. 2008), as well as to produce 
NAGase involved in chitin degradation (Miller et al. 
1998; Burke et al. 2011). 

Our results suggest that the complex interactions 
between the biotic and abiotic site factors may further 
propagate in affecting the soil microbial community func- 
tion involved in SOC transformation and/or turnover 
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through modification of the bacterial and saprophytic 
fungal community structure, shifts in the allocation of 
plant-derived organic carbon to microbial biomass, and 
regulation of microbial activities. Moreover, factors affect- 
ing soil microbial community structure may also impose 
direct and/or indirect impacts on the microbial function, 
enhancing or weakening the overall controls of soil 
microbial community function by individual site factors. 

We used the specific activity of the soil extracellular 
enzymes (i.e., enzymatic per unit SOC) to assess the 
dynamics of SOC transformation and/or turnover in path 
analysis. Soil bacteria community was identified to 
strongly regulate SOC storage although increasing the 
microbial C acquisition activity. Saprophytic fungi com- 
munity contributes to SOC turnover by increasing the 
microbial C oxidative activity. Extracellular enzyme activ- 
ity in soils has been studied for more than a century with 
a goal of understanding the biochemistry of decomposi- 
tion and nutrient cycling. Extracellular enzymes are 
widely produced by different groups of soil microorgan- 
isms, and some of them can be used as indicators of the 
presence or activity of specific microbial taxa (Baldrian 
2009). They provide useful insights into the mechanisms 
of microbial sensitivity to environment change (Cusack 
et al. 2011). The most widely assayed enzymes are those 
involved in the degradation of cellulose and lignin, the 
most abundant components of plant litter (Cusack et al. 
2011). Other commonly measured enzymes hydrolyze 
proteins, chitin, and peptidoglycan, which are the princi- 
pal reservoirs of organic N (Caldwell 2005). However, 
none of the soil enzyme assay techniques is able to com- 
pletely reflect the in situ conditions and to measure the 
actual rates of enzymatic reactions (Toberman et al. 
2008). Quantitative PCR of gene sequences or mRNA 
transcripts was recently used for quantification of some 
enzyme-encoding genes, including phenoloxidase-laccase 
(Hassett et al. 2009; Lauber et al. 2009). However, the 
connection between transcript or gene copy level and 
enzyme activity is indirect and is biased by genes encod- 
ing inactive molecules (Baldrian 2009). 

In summary, our study reveals complex controls of soil 
carbon dynamics by various above- and below-ground 
processes and their interactive effects. Soil microbial com- 
munity composition is strongly influenced by water, tem- 
perature, SOC, fine root mass, clay content, and C/N 
ratio in soils. Soil bacterial communities are strongly 
linked with the extracellular enzymes involved in carbon 
transformation, whereas saprophytic fungi are associated 
with activities of extracellular enzymes driving carbon 
oxidation. There exist complex interactions and linkage 
among plant traits, microenvironment, and soil physio- 
chemical properties in affecting SOC via microbial regula- 
tions in temperate forests. 
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